A method for determining event-by-event elliptic flow fluctuations 
based on the first-order event plane in heavy-ion collisions 
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A new method is presented for determining event-by-event fluctuations of elliptic flow, V2, us- 
ing first-order event planes. By studying the event-by-event distributions of V2 observables and 
first-order event-plane observables, average flow (t^) and event-by-event flow fluctuations can be 
separately determined, making appropriate allowance for the effects of finite multiplicity and non- 
flow. The method has been tested with Monte Carlo simulations. The connection between flow 
fluctuations and fluctuations of the initial-state participant eccentricity is discussed. 
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I. INTRODUCTION 



In heavy-ion collisions, the azimuthal distributions of 
emitted particles can be decomposed with a Fourier ex- 
pansion |l|: 



dN 

dip 



— {1 + 2v n cosn(^ - * RP )}, 



(1) 



where tp denotes the azimuthal angle of the particle and 
is the reaction plane azimuth (defined by the impact 
parameter vector). The Fourier coefficients, 



v„ = (cosn((p - *rp)) , 



(2) 



are referred to as anisotropic flow of the n th harmonic. 
The second harmonic, elliptic flow, carries information 
on the early stage of heavy- ion collisions, and has been 
extensively studied. Event-by-event flow fluctuations 
HI H, EL IE IE Q are °f considerable interest as they must 
be sensitive to the physics of the very early stages of 
the collision, and because any fluctuation observable has 
potential relevance for phase transition phenomena. Un- 
derstanding flow fluctuations would also greatly improve 
anisotropic flow measurements that are currently domi- 
nated by systematic uncertainties in which flow fluctua- 
tions play a crucial role. At fixed ccntrality, initial-state 
fluctuations in the spatial anisotropy of the participant 
zone will cause flow fluctuations. In addition to this 
inevitable source of fluctuations, there might be addi- 
tional fluctuation contributions that could offer unique 
insights into dynamical details of the collision process 
at very early times (1 fm/c and earlier) [f|. The large 
observed elliptic flow at the Relativistic Heavy Ion Col- 
lider (RHIC) points to a very short thermal equilibration 
time in the framework of hydrodynamic models. This 
puzzling feature calls for further investigation and could 
have alternative explanations or might be explained by 
exotic phenomena [1, IE HE [IH- Improved methods to 
experimentally determine flow fluctuations would be an 
important step towards addressing some or all of the open 
issues discussed above. 



Fluctuations in the shape of the initial participant re- 
gion, and in particular, in the orientation of the region's 
principal axes relative to the direction of the impact pa- 
rameter, lead to a non-trivial picture of anisotropic flow. 
In this picture, the apparent flow at mid-rapidity might 
be different in direction and magnitude from the real flow 
as measured with respect to the reaction plane. The 
method proposed in this paper is sensitive to such a dif- 
ference, and allows fluctuations in the orientation of the 
principal axes of the participant region to be measured. 

RHIC data hold much promise for the purpose of un- 
derstanding elliptic flow fluctuations, since v-i is large, 
while the statistical noise arising from finite multiplicity, 
that tends to obscure the dynamical fluctuations of inter- 
est, is smaller than at lower energies. Most flow analyses 
at RHIC to date have relied on the second-order event 
plane, whereas in the present study, a case is presented 
for utilizing the first-order event plane to determine the 
mean elliptic flow, and to isolate the sought-after dynam- 
ical fluctuations about that mean. In RHIC experiments, 
first-order event planes can be obtained, for example, via 
the ZDC-SMD (Zero Degree Calorimeter Shower Maxi- 
mum Detector) [H or the Forward TPC [U[ of the STAR 
detector. In the scenario envisaged here, the fluctuat- 
ing anisotropics are based on measurements near mid- 
rapidity, while the first-order event plane determination 
utilizes detectors that are far removed in rapidity. Con- 
sequently, non-flow effects (defined as azimuthal correla- 
tions that may contribute to v n measurements, but which 
are unrelated to the reaction plane orientation, or more 
generally, are unrelated to the initial geometry of the 
system) are believed to be negligible using this method 



II. TECHNIQUE 

With two independent first-order event planes ip a and 
ipb, elliptic flow can be determined with the help of the 
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relations 

V° hs = (C0S(2<^ - - *M) 

= (cos(2^ - 2* R p)) (cos(2* RP -ip a - Tpb)) 
= v 2 (cos(V> - ip b )) , (3) 

where the last factor, {cos(ip a — *Pb)) = (cos(ip a — 
^Rp))(cos(tpb — ^rp)) is theproduct of the two first- 
order event plane resolutions The above is based on 
the assumptions that the two event planes are indepen- 



dent, and that the distributions of ip a and ipb with respect 
to the true reaction plane are symmetric. 

We introduce two event-by-event quantities 

c 2 , k = (cos2(^-* RP;fe )) (4) 
S2-,k = (sin2(<£ - *Rp ;fc )} , (5) 

where index k denotes the fc th event and the average is 
taken over all particles in that event. Using the equality 



C0S(2<P - 1p a;k - 1p b;k ) = C0S[2(<^ - *RP ;fc ) - (tpa;k - *RP;fc) - {lpb;k ~ *RP;fc)] 

= cos2((p - f R p ; fe) cos(A^a;fc + Aipb-k) + sin 2(</j - ^rp^) sin(A^ a; fc + Aip b;k ), (6) 

where Aip a;k = ip a;k - *Rp ;fc and Aip b . k = ipb-k ~ ^rp^ (and a similar expression for sm(2(p - ip a . k - ip b - k )), and 
averaging over particles in the k th event, one finds 

C 2?k = c 2:fc C0s(A^ a; fc + Alp b - k ) + S 2 -k Sm(A?/'a;fc + Alp b - k ) (7) 

s 2 h k = s 2-,k cos(A-0a ; fc + Aip b ;k) - c 2 - k sm(Atp a . k + Aip b . k ) . (8) 



r 



We further assume that the distribution of the first-order 
event planes {Aip a . k + Aip b . k } is independent of distribu- 
tions {c 2 - k } and {s 2 - k } . This assumption is usually not 
valid for second-order event planes, and this illustrates 
one of the advantages of the present approach. Assuming 
that {Aip a . k } and {Aip b . k } are both symmetric around 
zero and independent, the distributions {Aip a -k + Aipb;k] 
and {Atp a . k — Aipb;k = Atp a b-k} are identical. We discuss 
the validity of all of the above assumptions in more detail 
at the end of Section [TTTT 

From Eqs. ([7]) and one can calculate the mean (now 
averaged over all events) and the mean square of and 



„obs 
b 2:k 



i c 2?k) = ( c 2;fc) (COS Alp ab[k ) = V 2 (COS A%p ab - k ) , (9) 



(^=0, ( 10 ) 



The fluctuations, <J 2 {c 2 -, k } = (c 2 . k ) — (c 2 - k ) 2 and 

& 2 {s 2 ;k} = (s 2 . k ), each have several contributions: dy- 
namical flow fluctuations, non-flow, and a statistical part 
that is related to finite event multiplicity, 



cr 2 {c 2;fc } = crl yn {c 2 -, k } + - + 



1 + 1-4 - 2(c 2;fe } 2 - S 



<J 2 {s 2 :k} = erL n {s 2 ;fe} + - 



2M 



l-v A 



2M 



(15) 



(16) 



where M denotes multiplicity, and S stands for the non- 
flow contribution. The v A term in the above equations 
arises from setting n = 2 in the equalities 

(cos 2 n((p - *rp)) - 0.5(1 + wan), (17) 
(sin 2 n(ip - *rp)> = 0.5(1 - v 2n ) . (18) 



((<%m = (<4 ;k ) (cos 2 A^ h;fc ) + (s 2 2 . k ) (sin 2 A^ ab . k ) 



(11) 

((stk) 2 ) = (4;k) (sin 2 A4^ k ) + {s% k ) (cos 2 A^ ab , k ) . 

(12) 

Conventionally, (cos Atp ab . k ) in Eq. ([9]) is regarded as a 
correction for the event plane resolution. From the above 
equations, one finds relations 

(ci k ) + (sl k ) = {(c£ k s ) 2 ) + {(s?£) 2 ), (is) 



/ 2 \ / 2 \ (( C 2:A;) ) (( S 2;fe ) ) 

( c 2:k) ~ ( S 2:k) = > 0A , C ■ (I 4 ) 

(cos2A^ ab;fc ) 



Note that V4 is usually negligible compared to 1. Also, 
the terms inversely proportional to multiplicity can be 
experimentally measured by studying the dependence of 
c 2 {c2;fe} and <r 2 {s 2 - k } on the multiplicity of particles 
used in the event. Either measuring these terms, or just 
neglecting the difference in ~ 1/M terms for a 2 {c 2 . }k } 
and <7 2 {s2 ; fc} and evaluating the difference cr 2 {c 2 - k } — 
c 2 {s2;fc}, one gets access to the difference cr^ yn {c 2 - k } — 
cr 2 yn {s2;fc}. The latter is directly related to flow fluctua- 
tions, but as shown below, it also depends on fluctuations 
in the orientation of the principal axes of the partici- 
pant region with respect to the direction of the impact 
parameter. We return to this question after discussing 
eccentricity fluctuations in Section III. 
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Relations established above have been tested in a 
Monte Carlo simulation. In each event, the azimuthal 
angle of each particle at mid-rapidity has been assigned 
randomly according to the distribution of Eq. JT]). In 
this Monte Carlo simulation, only vi is non-zero, and 
non-flow effects have not been implemented. From event 
to event, the vi value fluctuates according to a Gaussian 
distribution. The first-order event plane follows a typical 
event plane distribution described in Ref [l[. In the first 
set of simulations, the input corresponds to 5% mean el- 
liptic flow and 3% dynamical flow fluctuations. We set 
each of the first-order event plane resolutions to be 20%, 
which corresponds to (cos Aip a i,- k ) = 4%. Five multiplic- 
ities are tested: M = 25, 50, 100, 200 and 400. In each 
case, ten million events were generated. Table H] lists the 
output results. Over a broad range of multiplicities, the 



Multiplicity 




0-dyn{c 2 ;fc}(%) 


25 


4.90 ±0.12 


3.16 ±0.18 


50 


4.97 ±0.09 


3.00 ±0.14 


100 


4.90 ±0.07 


3.14 ±0.11 


200 


4.96 ±0.06 


3.03 ±0.10 


400 


5.05 ±0.05 


2.94 ±0.09 



TABLE I: Results for the extracted mean elliptic flow for an 
input of 5%, and its dynamical fluctuation for an input of 3%, 
based on simulated events with different multiplicities. The 
errors are statistical. 

output results well agree with the input values within the 
statistical errors of the simulation. 

We have also explored the robustness of the method 
with variations in the input values of (c2-.k) and 
Cdyn{c2;fe}- In this second group of tests, the multiplicity 
was fixed at M = 100, and one million events were gen- 
erated in each case. The reconstructed mean V2 and the 
extracted dynamical fluctuations, Od yn {c2 ; fc}, are shown 
in Fig. [T] and they are found to be consistent with the 
input values. 



III. ECCENTRICITY FLUCTUATIONS 

In heavy-ion collisions, due to the finite number of par- 
ticipants, the center of the overlap zone can be shifted 
and the orientation of the principal axes of the interac- 
tion zone can be rotated with respect to the conventional 
coordinate system with the x axis pointing along the im- 
pact parameter, as illustrated in Fig.O As a result, the 
final particle distribution is symmetric about the x' axis, 
instead of the x axis. The effect of the shift is negligi- 
ble compared to the effect of the rotation, and therefore 
we concentrate on the latter. The conclusion that the 
shift can be neglected is based on detailed Monte Carlo 
calculations using a Glauber model where the geometry 
of the interaction zone is defined by the position of the 
participating nucleons. 




<c 2 / n, = 4% 




<C 2:k > (i V% 



1 1.5 2 2.5 3 3.5 1 1.5 2 2.5 3 3.5 
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FIG. 1: (color online) The reconstructed mean V2 (stars) 
and its dynamical fluctuation (circles) as functions of input 
Odyn{c2 ; fc} for Monte Carlo events with multiplicity of 100. 
The errors are statistical. Panels are shown for input (c2 ; fc) 
ranging from 2% to 5%. The horizontal lines represent the 
input mean V2, and the diagonal lines represent the input 
dynamical fluctuation of V2- 



We define the angle between the x' and x directions to 
be Atp e . k for the fc th event. 

The standard definition of eccentricity [i| for the k th 
event is 



(Vi ~ x l)k 
{Vi + 



(19) 



where xi and denote the i-th participant coordinates 
and the average is taken over all participants in the event. 
The coordinates (x,y) and (x',y') are linked by the ro- 
tation through Aip £ -k, 

x k = x' k cos Aip £ . k - y' k sin Aip e . k 

y k = x' k sin Aip E . k + y' k cos Aip E . k , (20) 

which leads to the relation 



s k = 4 cos2AV> £;fc . 



(21) 



Note that from the definition of the rotated frame, 
( x 'i Hi) = 0- The elliptic anisotropy in a given event is 
developed in the x',y' plane such that {s' 2 . k ) = and 
<7,}yn{s2 fc} = 0- Also, similar to the relation between 
eccentricities, one finds 



(C 2 ;fe) = (4;fc) (cOs2AV'£;fe) , 

which follows directly from 



C2;fe 
S2:k 



c' 2 . k cos2Aip e . k 
c' 2k sm2Aip e . k 



s' 2 . k sin 2 Aip £ , k 
s' 2 . k cos2Aip Elk 



(22) 

(23) 
(24) 
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FIG. 2: (color online) Schematic representation, in the plane 
transverse to the beam (z) direction, of a collision between 
two identical nuclei. The x- and y-axes are drawn as per the 
standard convention. The solid circles illustrate a possible 
configuration of the participant nucleons. Due to the fluctu- 
ation of this event, the overlap zone is shifted and tilted with 
respect to the (x, y) frame, x' and y' are the principal axes 
of inertia of the solid circles. 

Therefore (c2-k) is always less than or equal to (c' 2 . fc ). 
Note that 

4fc + 4fe = 4? fc + «affc . (25) 

which reflects the fact that such a combination is inde- 
pendent of the plane in which it is calculated, as it de- 
pends only on the particle pair angle differences in the 
event. As mentioned in the previous section, this sum 
contains the flow fluctuation contribution (including the 
fluctuation in the orientation of the principal axes of the 
participant zone) as well as the non-flow contribution. 
The non-statistical part of this sum (everything except 
~ 1/M terms) corresponds exactly to u|{2} — ellip- 
tic flow measured with two-particle correlations at mid- 
rapidity. To remove non-flow contributions, one should 
consider the difference 

<4* " 4k> = ( c 2% - *2? fc > <cos4AV> £;fc >. (26) 
The non-statistical part of this difference, 

(4;k - S 2;fc)dyn = (c' 2 %)dyn (cOs4A^ e;fc ), (27) 

provides an important relation between fluctuations mea- 
sured with respect to the first-order reaction plane, 
( c 2-k ~ s 2-fe)dyn, flow fluctuations measured at mid- 
rapidity (which includes effects of the fluctuations in the 
geometry of the participant zone), and the distribution 
in A?/> e ;fc. 

In the picture described above, when the flow fluctu- 
ations are driven by fluctuations in the participant ec- 
centricity, it is not obvious that the two first-order event 



planes defined by spectators from the two nuclei are inde- 
pendent, nor is it obvious that the first-order event plane 
is independent of the second-order event plane defined by 
the participants. Indeed, the positions of the spectators 
are somewhat correlated with the positions of the partic- 
ipants, but as we found using the Monte Carlo Glauber 
model, this has a negligible effect on correlations of the 
event planes. In our study, we used the center of gravity 
of the spectator distribution with respect to the nuclear 
center to define the first-order event plane for each nu- 
cleus (ipi,a and '01, b) i an d the second-order event plane 
(1P2) was defined by the minor axis of the participant 
zone. Using the center of the collision instead of the cen- 
ter of gravity of the spectator distribution would make 
the correlation effects even smaller. 

We find that for most centralities, the correlations 

(COS(V>1, " *RP)) (C0S(V>1,& - *Rp)) 

and 

(cos(y>i, a + -0i, b - 2tp 2 )} 

(COS(V>1, - *RP)> (C0S(V>1,& - 4-Rp)) <C0s(2l/>2 ~ 2* RP )) 

-1 (29) 

are at the sub-percent level, with a maximum of about 
a few percent for the 5% most central collisions. Be- 
sides the event plane correlations due to the correlated 
positions of spectators and participating nucleons, mo- 
mentum conservation also deserves consideration. We 
refer here to an experimental study [Til ] which concluded 
that there is negligible momentum-conservation correla- 
tion between the event plane based on spectators from 
each nucleus separately and the orientation of directed 
flow close to mid-rapidity. 

IV. DISCUSSION AND SUMMARY 

It is useful at this stage to consider the notation 
(c' 2 . k ) = (v' 2 ) and (c2 2 fc )dyn = W2 2 ) ■, where v' 2 is the "ap- 
parent" flow at mid-rapidity — elliptic event anisotropy 
measured with respect to the principal axes of the partici- 
pant zone. Much progress towards a direct measurement 
of flow fluctuations at mid-rapidity has been reported 
recently [f| Then the above equations can provide 
important information on fluctuations in the orientation 
of the principal axes of the participant zone. 

Note that the orientation of the participant zone (and, 
consequently, the apparent anisotropic flow) can depend 
on the rapidity of the particles under study. Then, in 
principle, one can study the correlations in the orienta- 
tion of anisotropic flow as function of particle rapidity. 
Such information will be very valuable for the reconstruc- 
tion of the initial conditions in heavy-ion collisions. 

In summary, various suggestions in the literature point 
to dynamical event-by-event fluctuations in elliptic flow 
as being of great interest in the realm of RHIC physics, 
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and such fluctuations are argued to be especially relevant 
for understanding collision dynamics at the earliest times. 
On a more technical level in experimental methodology, 
the magnitude of flow fluctuations associated with cur- 
rent elliptic flow measurements is not understood, and 
this uncertainty affects the overall systematic error on 
these measurements. Prompted by the above consider- 
ations, this work presents a new method for experimen- 
tal analysis of elliptic flow in a scenario where the first- 
order event plane can be resolved. The method allows 
the extraction of mean V2 and its dynamical cvent-by- 
event fluctuations, and good immunity to both statistical 
fluctuations and non-flow effects can be expected. Sim- 
ulations have been presented that validate the method 



under a range of conditions similar to those observed in 
RHIC data, ft has been shown that measurements of flow 
fluctuations using the first-order event plane, accompa- 
nied by measurements of apparent flow fluctuations at 
mid-rapidity, can also provide important information on 
the fluctuations of the participant zone. 
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